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I. INTRODUCTION 

Light-cone Hamiltonian diagonalization methods offer a number of attractive advantages 
for solving nonperturbative problems in quantum field theory, such as a physical Minkowski 
space description, boost invariance of the bound-state wave functions, no requirement for 
fermion doubling, and a consistent Fock-state expansion well matched to physical problems. 
In the discretized light-cone quantization (DLCQ) method, the light-cone Hamiltonian H^c 
of a quantum field theory is diagonalized on a discrete Fock basis defined by assuming peri- 
odic boundary conditions in the light-cone coordinates |],@]. The eigenvalues of H L q give the 
mass spectrum of the theory, and the respective eigenfunctions projected on the free Fock 
basis provide the frame-independent light-cone wave functions needed for phenomenology || 
including the amplitudes needed to compute exclusive B decays ||-|5| , deeply virtual Comp- 
ton scattering |7|||, and other hard exclusive processes 0. The DLCQ method has been 
successfully used to solve a large number of one-space and one-time theories [Q], including 
super symmetric gauge theories |I0| . It also has found application in analyzing confinement 
mechanisms [IT], string theory |12|| , and M-theory [|13[1 . 



The application of DLCQ to physical, (3+l)-dimensional space-time quantum field the- 
ories is computationally challenging because of the rapid growth of the number of degrees of 
freedom as the size of the Fock representation grows. A promising alternative is the trans- 



verse lattice method [I i] which combines light-cone methods in the longitudinal light-cone 



direction with a space-time lattice for the transverse dimensions. Recently Dalley |I1| and 



Burkhart and Seal |I(J have extended the transverse lattice method to estimate the shape of 
the valence light-cone wave function of a pion, a key input to much hadron phenomenology. 
Burkhart and Seal have also given an explicit calculation of the Isgur-Wise function for 
semileptonic B decays fL7| . 

Another major difficulty in applying DLCQ to quantum field theory in 3+1 dimensions 
is the implementation of a nonperturbative renormalization method. Most methods of reg- 
ulating nonperturbative calculations in the light-cone representation, such as momentum 
cutoffs, do not allow a correct renormalization even of perturbative calculations. The prob- 
lem can be traced to the fact that any momentum cutoff violates Lorentz invariance as well 



as gauge invariance ||18|| . Since dimensional regularization is not available in DLCQ, one 
needs to introduce new fields or degrees of freedom to render the ultraviolet behavior of 
the theory finite. One intriguing possibility is to analyze ultraviolet-finite supersymmetric 
theories and then introduce breaking of the theory. The heavy supersymmetric partners 
then regulate the ordinary sector of the theory in a manner analogous to Pauli-Villars (PV) 



regulation |T9| . String theory also provides mechanisms for regulating quantum field theory 
at short distances which are equivalent to an infinite spectrum of PV particles [f2D[. The 
introduction of PV fields can thus regulate a theory covariantly, after which the discretized 
momentum grid of DLCQ acts only as a numerical tool in the manner of performing a 
numerical integral. 

In our previous work [TB,2T| we have shown that a model field theory in 3+1 dimensions 
can be solved by combining DLCQ with PV regulation of the ultraviolet regime. In our 
first application [|l8j, a model theory was constructed to have an exact analytic solution by 
which the DLCQ results could be checked, for both accuracy and rapidity of convergence. 
The model was regulated in the ultraviolet by a single PV boson, which was included in 



the DLCQ Fock basis in the same way as the "physical" particles of the theory. We then 
extended this approach to a more realistic model which mimics many features of a full 
quantum field theory |21|| . Unlike the analytic model which contained a static source, the 
light-cone energies of the particles in this model have the correct longitudinal and transverse 
momentum dependence. 

An important question is whether the generalized PV method with a finite number of 
fields can regulate a field theory at all orders. Paston and Franke E2] have studied the 



relation between perturbation theory in the light-cone representation and standard Feyn- 
man perturbation theory, and they have developed general rules for testing regularization 
procedures. For full Yukawa theory, Paston, Franke and Prokhvatilov |23|] have shown that 
one PV boson and two PV fermions can regulate the theory in such a way as to allow a 
correct perturbative renormalization. 

In this paper we shall apply generalized PV regularization and discrete light-cone quan- 
tization to the nonperturbative solution of (3+l)-dimensional Yukawa theory in a single- 
fermion truncation. We allow any number of bosons in the Yukawa theory but only one 
fermion in the Fock representation; fermion pair terms and any other terms that involve 
anti-fermions are neglected. We shall thus consider a field-theoretic model where one par- 
ticle, which we take to be a fermion of mass M, acts as a dynamical source and sink for 
bosons of mass \x. In addition, three heavy PV scalars, including two with negative norm, 
will be used to regulate the theory so that the chiral properties of the renormalized theory 
are maintained, at least to one loop in perturbation theory. In particular, the mass of the 
renormalized fermion constituent vanishes as its bare mass vanishes. A distinct advantage of 
our approach is that the counterterms are generated automatically by the PV particles and 
their negative-metric couplings. We emphasize that the PV fields are added from the be- 
ginning, at the Lagrangian level, to facilitate our nonperturbative calculations, rather than 
being invoked for subtractions at a diagrammatic level, which would limit the implemen- 
tation to perturbation theory. Note that the PV fields enter singly, doubly or triply at all 
three-point vertices. This contrasts with use of generalized Pauli-Villars spectral integrals 
over mass to regulate divergences. 

The extra degrees of freedom of the P V scalars and their negative-metric couplings intro- 
duce new computational challenges. However, the matrix eigenvalue problem can be solved 
for the lowest-mass state by the use of a new, indefinite-metric Lanczos algorithm which we 
describe in an Appendix. We also calculate the light-cone wave function of each Fock-sector 
component and the values for various physical quantities, such as average multiplicities and 
average momenta of constituents, bosonic and fermionic structure functions, and a form 
factor slope. We also verify that with our choice of PV conditions, the DLCQ calculations 
of the nonperturbative theory at weak coupling coincide with the covariant perturbation 
theory through one loop, although numerical resolution does start to become a problem 
even on a supercomputer. 

In our convention, we define light-cone coordinates pl| by 

x ± = x° + x 3 , x±— (x 1 ,x 2 ) . (1.1) 

The time coordinate is taken to be x + . The dot product of two four- vectors is 

p ■ x = ~(p + x~ + p~x + ) — p± ■ x± . (1.2) 



Thus the momentum component conjugate to x is p + , and the light-cone energy is p . We 
use underscores to identify light-cone three-vectors, such as 



P= (P + ,P±)- 



;i.3) 



For additional details, see Appendix A of Ref. ||T8[ or a review paper 



The following is an outline of the remainder of the paper. In Sec. [II] we discuss the 
regularization and renormalization of the Yukawa Hamiltonian. Our numerical methods 
and results are presented in Sec. fTJ. Section [TV] contains our conclusions and plans for 
future work. Details of the numerical diagonalization method are given in the Appendix. 



II. YUKAWA THEORY 



A. Light-cone Hamiltonian 



We write the Lagrangian for Yukawa theory as 



c = \{d^f-\^ + l - 



Th^dp. - (<WV ^ - Mipip - gcfripip . 



(2.1) 



The corresponding light-cone Hamiltonian has been given by McCartor and Robertson [25 



Here we consider a single-fermion truncation and therefore neglect pair terms and any other 
terms that involve anti-fermions. We also neglect longitudinal zero modes. To regulate the 
theory, we include three PV bosons. The resulting Hamiltonian is 
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+ h.c. 



+&L»^A««L%m (^ _ ^^ + ^ + m)/ 



where M is the fermion mass, \x = /io is the physical boson mass, //j and (—1)* are the mass 
and norm of the i-th PV boson, and e<i s = — js(2s,i). The nonzero commutators are 



a im,Q'j m i — \ J- J VijVm,m' > i "n,s; "n',s' J — "n,n'"s,s' ■ [.*•**) 



The different boson couplings are denoted by ^g, where £ — 1 corresponds to the physical 
boson; the other £j are chosen to arrange cancellations discussed below. Fermion self-induced 
inertia terms cancel in a sum over bosons and therefore do not appear. The bare parameters 
are the coupling g and the mass shift 5M 2 . 



The number of PV flavors is determined by the cancellations needed to regulate the 
theory and restore chiral invariance in the M = limit [^,|l^]. The discrete chiral symmetry 
ip — > i^ip, 4> — > —0 should itself guarantee a zero self-mass in this limit, if not for the fact 
that the symmetry is typically broken by a cutoff. The one- loop self-energy of the fermion 

is m 
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with A 2 a cutoff such that A 2 ^> /i 2 ^> M 2 . In order that the self-energy be finite and 
proportional to M 2 , the relative coupling strengths £, must satisfy the constraints 



1 = 1 1=1 1 = 1 



0. 



(2.5) 



In addition, the norm of the i-th PV field must be chosen as (—1)*. 

The third constraint in (|2.5|) is peculiar to the one-loop calculation. For higher-order 
or nonperturbative calculations, it must be replaced by a more general renormalization 
condition. However, to simplify the numerical work, we use the one-loop constraint and 
then check for failure of cancellation in the zero- mass limit. 

As it stands, H^q contains infrared divergences associated with the instantaneous fermion 
interaction, which is singular at the point where the instantaneous fermion has zero longi- 
tudinal momentum. The divergences are partly cancelled by crossed-boson contributions. 
To cancel the remainder we need to add an effective interaction, modeled on the missing 
fermion Z graph. The effective interaction is constructed from the pair creation and annihi- 



lation terms in the Yukawa light-cone energy operator [25[ 
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The denominator for the intermediate state is 

M 2 _ _ M 2 + p'l M 2 + (q' ± - p ± ) 2 _M 2 +p 2 ± 
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(2.7) 



To guarantee the cancellation of the singularity in the numerical calculation, the instanta- 
neous interaction is kept only if the corresponding crossed-boson graph is permitted by the 
numerical cutoffs. 

The Fock-state expansion for the single-fermion eigenstate of the Hamiltonian is 



.,=v«p^ t i^nim^Y. 



no,ni,n2,ri3=0 
ntot 



y/lQn' 6 ^ 



'lGiT 3 qj 

ntot 1 "tot 



(2.8) 



"tot 



where no is the number of physical bosons, rii the number of PV bosons of flavor i, n tot = 
Si=o n i' an d h i s the flavor of the j-th constituent boson. It solves the eigenvalue problem 
Hhc^a — M 2 $ CT . The normalization is 



&} ■ $ ff = 1Qti 3 P + 5(F - P) . 



(2.9) 



B. Renormalization conditions 

Mass renormalization is carried out by rearranging the eigenvalue problem into one for 
5M 2 at fixed M 2 
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Here x = p + / P + is the fermion momentum fraction, /C is shorthand for the interaction 
kernel, and the (f) = <p/ 'y/x are new wave functions. 

The coupling is fixed by setting a value for the expectation value (:</> 2 (0):) = $J. :0 2 (O): Q a 
for the boson field operator (p. This quantity is useful because it can be computed fairly 
efficiently in a sum similar to a normalization sum 
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with (— l)( n ») being the norm of the state with boson flavor partitioning (n. 



III. RESULTS 



A. Numerical methods 

The principal tools for the solution of the Hamiltonian eigenvalue problem are DLCQ || 
and the Lanczos diagonalization algorithm [27|. The first converts the problem to a matrix 
form, and the second quickly generates one or more eigenvectors. The use of negatively- 
normed states makes the ordinary Lanczos algorithm inapplicable; however, an efficient 



6 



generalization has been developed for this situation. It is discussed in the Appendix. The 
constraint for the coupling renormalization is solved iteratively. 

The discretization is based on the standard DLCQ approach where longitudinal and 
transverse momenta are assigned discrete values q + = rmr/L and q ± = n±7r/Lj_, with L and 
L± chosen length scales. Momentum conservation requires the individual m to sum to a fixed 
constant K, where P + = Ktc/L is the total longitudinal momentum. The integer K is called 
the harmonic resolution [Q], because longitudinal momentum fractions y = q + /P + = m/K 
are resolved to order 1/K. The positivity of longitudinal momenta forces a natural cutoff 
such that m < K. Also, the eigenvalue problem for M 2 is independent of L; the length scale 
cancels between P + and P~ . 

The boundary conditions in the longitudinal direction are chosen to be periodic for the 
boson fields but antiperiodic for the fermion field. This means that the integers m are even 
for bosons, and the corresponding fermion momentum index is odd. The harmonic resolution 
K is then also odd for the single-fermion state considered here. 

The transverse direction requires an explicit cutoff A 2 , which we impose on individual 
light-cone energies 

(fi 2 + m 2 ± 7r 2 /L 2 ± )/m< A 2 /K. (3.1) 

The total transverse momentum is taken to be zero. The integers n x and n y are limited to 
a range [— iVj_, N±\, which, along with the cutoff, determines the transverse scale L±. 

Given this discretization, the eigenvalue problem is converted to a matrix problem by 
a trapezoidal approximation to any momentum integrals. We have found useful modifi- 
cations jn],[n|] which include non-constant weighting factors near the integral boundaries. 
These weights are adjusted to compensate for the DLCQ grid being incommensurate with 
the boundaries. In the present calculation, these weights are kept only in the two-body 
sector where maximal symmetry can be maintained. For higher Fock sectors, sensitivity 
to cancellations is of greater concern than boundary effects, and the weights disrupt the 
cancellations. 

Unlike an ordinary eigenvalue problem, the presence of negatively normed states allows 
unphysical states to be lower in mass than the physical one-fermion state. Criteria must 
be employed to select the correct state in a numerical calculation. We used the following: 
a positive norm, a real eigenvalue, a minimum number of nodes (preferably zero) in the 
parallel-helicity boson-fermion wave function, and the largest bare-fermion probability be- 
tween and 1. Each of these characteristics can be computed without constructing the 
full eigenvector, provided one saves a few components of each Lanczos vector g- (see the 
Appendix). 

As a check on the calculation, we took advantage of an exact solution that exists for 
the unphysical situation of equal- mass PV bosons Iffijjl . In a particular (null) basis the 
matrix representation is purely triangular. Each wave function of the dressed fermion is 
then directly computable in a finite number of steps. 

For comparison, we also solved the problem using Brillouin-Wigner perturbation theory, 
for which 

_ / i~^io)(o|6g y , 



The integrals were computed numerically with the same discretization as the nonpertur- 
bative DLCQ calculation. The main effort in the perturbative expansion is then matrix 
multiplication, just as for the Lanczos algorithm. 



B. Computed quantities 

Various quantities can be computed from the wave functions (fft 1 *' for the different Fock 
sectors. We compute the slope of the no-flip form factor of the fermion, structure functions 
for bosons and the fermion, average momenta, average multiplicities, and a quantity sensitive 
to boson correlations. The form factor slope F'(0) is given by 
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The physical boson structure function for bare helicity s is defined as 
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The normalization is such that the integral yields the average multiplicity (ub)- We separate 
two pieces, for parallel and antiparallel bare helicity 



(3.5) 
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and have 

(n B ) = (n B+ ) + (n B -) ■ 
The average boson momentum is treated analogously, with 

(Vb s ) = / yf Bs (y)dy and (y B ) = (y B+ ) + (y B _) . 
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As a measure of the correlations in the multiple-boson Fock sectors, we compute ((3/13/2) — 

(y) 2 )n>2/(y) 2 where 
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and (y) n >2 is the same as (y) except that only states with two or more bosons are included. 
Calculations at low resolutions tend to have difficulty for stronger coupling. This can be 
seen already at third order in perturbation theory, where loop integrals are poorly approxi- 
mated and subtractions between loops with different boson flavors magnify the errors. Fully 
averaged quantities such as (tib) are less affected by this, but the structure function Jb can 
be quite poorly represented. An example is given in Fig. [l]. 
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FIG. 1. The boson structure function fg at various low to moderate numerical resolutions, 
with M = fi, (:0 2 (O):) = 1.0, A 2 = 50^ 2 , //? = 10/i 2 , /j% = 20/x 2 , and /i§ = 30/i 2 . To compare with 
a structure function computed at higher resolutions, see Fig. |2|( 



Clearly one cure is to work at higher resolution. Limitations on computer storage then 
require truncation in the number of constituents. Luckily, even at the strongest coupling 
that we consider, states with large numbers of constituents are unimportant, and yet the 
results differ significantly from lowest-order perturbation theory. The relative importance of 
different numbers of constituents can be seen in Table |, where we list probabilities for the 
various Fock sector contributions to a calculation with moderate resolution. Truncation to 
a maximum of two bosons is seen to offer a very good approximation. For weaker couplings, 
lower resolution is sufficient. 

Given these considerations for resolution and truncation, we have done two sets of cal- 
culations with K between 11 and 29, or even 39. For K = 11 and 13 there is no explicit 
truncation and the maximum number of bosons is 5 and 6, respectively. For K — 15, 17, 
and 19, the maximum number of bosons used is 4, and for K > 21 the maximum used is 
2. Two different sets of cutoff and PV masses were considered: A 2 = 50/i 2 , /i 2 = 10/i 2 , 
fj% = 20/x 2 , and f4 = 30/i 2 ; and A 2 = 100/i 2 , /i 2 = 20/i 2 , \i\ = 40/i 2 , and fj% = 60/x 2 . The 
transverse resolution A^j_ was at least 4 and was increased beyond that to the extent allowed 
by the available storage on a 16 GB node of an IBM SP. The four processors of the node 
were used in parallel. 

A sampling of explored parameter values can be seen in Tables [TT[ jVll| . For each choice 
of input parameter values, these tables present the results for the bare parameters of the 
Hamiltonian, g and 6M 2 , and for various expectation values. 



TABLE I. Fock sector probabilities f | 
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the number of physical bosons and rii the number of Pauli-Villars bosons of type i. The helicities a 
and s refer to the physical and bare fermion, respectively. The numerical and physical parameters 
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TABLE II. Bare parameters and observables. The input parameter values are M 2 
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0.097 


2868. 


19 


4 


0.100 


1.287 


0.054 


0.988 


0.091 


3324. 


19 


5 


0.100 


1.349 


0.056 


0.988 


0.096 


3375. 


21 


4 


0.100 


1.281 


0.055 


0.988 


0.088 


4075. 


21 


5 


0.100 


1.350 


0.057 


0.988 


0.095 


3257. 


21 


6 


0.100 


1.378 


0.058 


0.989 


0.099 


2788. 


21 


7 


0.100 


1.392 


0.058 


0.989 


0.101 


2977. 


21 


8 


0.100 


1.400 


0.059 


0.989 


0.104 


2920. 


29 


4 


0.100 


1.302 


0.062 


0.987 


0.091 


3572. 


29 


5 


0.100 


1.353 


0.063 


0.988 


0.096 


3441. 


29 


6 


0.100 


1.386 


0.065 


0.989 


0.099 


2920. 


29 


7 


0.100 


1.400 


0.066 


0.989 


0.103 


2983. 



TABLE III. Additional observables. The input parameter values are the same as for Table ||. 
The dashes represent values not extracted from preliminary, lower-resolution calculations. 



K 


N ± 


(:^ 2 (0):> 


( n B,a) 


(n s _ CT ) 


(ran) 


{VB,<t) 


(VB,-o) 


(vb) 


11 


4 


0.100 


- 


- 


0.021 


- 


- 


0.01191 


11 


5 


0.100 


- 


- 


0.021 


- 


- 


0.01209 


11 


6 


0.100 


- 


- 


0.021 


- 


- 


0.01213 


13 


4 


0.100 


- 


- 


0.021 


- 


- 


0.01170 


13 


5 


0.100 


- 


- 


0.022 


- 


- 


0.01219 


15 


4 


0.100 


- 


- 


0.021 


- 


- 


0.01198 


15 


5 


0.100 


- 


- 


0.022 


- 


- 


0.01225 


17 


4 


0.100 


- 


- 


0.022 


- 


- 


0.01220 


17 


5 


0.100 


- 


- 


0.022 


- 


- 


0.01234 


19 


4 


0.100 


- 


- 


0.022 


- 


- 


0.01231 


19 


5 


0.100 


- 


- 


0.022 


- 


- 


0.01238 


21 


4 


0.100 


0.0147 


0.0072 


0.0219 


0.00836 


0.00417 


0.01253 


21 


5 


0.100 


0.0133 


0.0086 


0.0218 


0.00735 


0.00506 


0.01241 


21 


6 


0.100 


0.0128 


0.0091 


0.0220 


0.00712 


0.00540 


0.01252 


21 


7 


0.100 


0.0125 


0.0094 


0.0219 


0.00695 


0.00556 


0.01251 


21 


8 


0.100 


0.0125 


0.0096 


0.0220 


0.00691 


0.00566 


0.01257 


29 


4 


0.100 


0.0145 


0.0076 


0.0221 


0.00808 


0.00444 


0.01252 


29 


5 


0.100 


0.0135 


0.0086 


0.0221 


0.00753 


0.00508 


0.01260 


29 


6 


0.100 


0.0129 


0.0092 


0.0221 


0.00711 


0.00546 


0.01257 


29 


7 


0.100 


0.0126 


0.0095 


0.0221 


0.00697 


0.00564 


0.01261 
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TABLE IV. Same as Table HI, but for stronger coupling. 



K 


N ± 


(:^(0):> 


g 


SAP/fj, 2 


l^ol 2 


-100[i 2 F'(0) 


({VlV2)-{v) 2 )n>2 

(vY 1 


11 


4 


1.000 


5.417 


0.863 


0.891 


1.340 


21.42 


11 


5 


1.001 


5.325 


0.699 


0.899 


1.073 


27.77 


11 


6 


1.000 


5.105 


0.819 


0.883 


2.849 


22.06 


13 


4 


1.000 


4.520 


0.362 


0.886 


0.929 


36.62 


13 


5 


0.998 


4.746 


0.362 


0.893 


1.244 


28.91 


15 


4 


1.000 


5.783 


1.161 


0.918 


0.457 


22.85 


15 


5 


1.000 


5.589 


1.124 


0.913 


0.641 


23.52 


17 


4 


1.000 


5.045 


0.619 


0.909 


0.683 


25.69 


17 


5 


0.998 


4.806 


0.562 


0.899 


1.032 


26.31 


19 


4 


1.000 


4.804 


0.576 


0.895 


0.853 


28.82 


19 


5 


1.000 


4.973 


0.597 


0.902 


0.868 


29.83 


21 


4 


1.000 


4.925 


0.615 


0.896 


0.774 


32.23 


21 


5 


1.000 


5.082 


0.616 


0.903 


0.876 


28.27 


21 


6 


1.000 


4.978 


0.619 


0.900 


1.266 


25.64 


21 


7 


0.999 


5.046 


0.612 


0.903 


1.140 


27.78 


21 


8 


1.000 


5.009 


0.622 


0.901 


1.248 


27.68 


29 


4 


1.000 


5.089 


0.713 


0.902 


0.659 


29.23 


29 


5 


1.000 


5.072 


0.695 


0.901 


0.837 


29.11 


29 


6 


1.000 


5.173 


0.708 


0.907 


0.867 


26.29 


29 


7 


0.999 


5.152 


0.712 


0.905 


0.933 


27.62 



TABLE V. Same as Table III, but for stronger coupling. 



K 


N ± 


<:</> 2 (0):> 


( n S,<7> 


i n B,-<j) 


<"B> 


(VB.a) 


{VB-a) 


(vb) 


11 


4 


1.000 


- 


- 


0.237 


- 


- 


0.1425 


11 


5 


1.001 


- 


- 


0.231 


- 


- 


0.1342 


11 


6 


1.000 


- 


- 


0.233 


- 


- 


0.1383 


13 


4 


1.000 


- 


- 


0.217 


- 


- 


0.1205 


13 


5 


0.998 


- 


- 


0.222 


- 


- 


0.1280 


15 


4 


1.000 


- 


- 


0.210 


- 


- 


0.1167 


15 


5 


1.000 


- 


- 


0.212 


- 


- 


0.1181 


17 


4 


1.000 


- 


- 


0.211 


- 


- 


0.1192 


17 


5 


0.998 


- 


- 


0.214 


- 


- 


0.1212 


19 


4 


1.000 


- 


- 


0.216 


- 


- 


0.1220 


19 


5 


1.000 


- 


- 


0.216 


- 


- 


0.1217 


21 


4 


1.000 


0.1179 


0.0994 


0.2173 


0.06625 


0.05649 


0.12275 


21 


5 


1.000 


0.1044 


0.1129 


0.2174 


0.05788 


0.06580 


0.12368 


21 


6 


1.000 


0.1031 


0.1154 


0.2185 


0.05770 


0.06759 


0.12529 


21 


7 


1.000 


0.1008 


0.1162 


0.2170 


0.05575 


0.06751 


0.12326 


21 


8 


1.000 


0.1016 


0.1153 


0.2169 


0.05608 


0.06655 


0.12263 


29 


4 


1.000 


0.1088 


0.1057 


0.2145 


0.05894 


0.06072 


0.11966 


29 


5 


1.000 


0.1048 


0.1122 


0.2170 


0.05738 


0.06452 


0.12190 


29 


6 


1.000 


0.0982 


0.1178 


0.2160 


0.05363 


0.06824 


0.12187 


29 


7 


1.000 


0.0985 


0.1180 


0.2165 


0.05360 


0.06803 


0.12163 
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TABLE VI. Same as Table 0, but for A 2 = 100/x 2 , ^ 2 = 20/i 2 , \i\ = 40/i 2 , and n\ = 60/j 



K 


N ± 


(:^(0):> 


g 


8M 2 /tj, 2 


l^ol 2 


-100^ 2 -F'(0) 


«!/l!/2>-<!/} 2 )„>2 

(y) 2 


11 


4 


0.989 


4.661 


1.220 


0.858 


1.195 


22.07 


11 


5 


0.998 


4.326 


0.991 


0.868 


1.114 


26.07 


11 


6 


0.904 


4.181 


1.063 


0.882 


1.442 


28.64 


13 


4 


1.002 


3.378 


0.487 


0.836 


1.122 


32.87 


13 


5 


1.000 


4.182 


0.537 


0.863 


1.245 


27.45 


15 


4 


1.000 


5.332 


1.849 


0.912 


0.337 


19.69 


15 


5 


1.000 


5.142 


1.804 


0.908 


0.458 


21.64 


17 


4 


0.998 


4.448 


1.131 


0.896 


0.448 


28.27 


17 


5 


1.000 


4.527 


1.094 


0.889 


0.758 


24.09 


19 


4 


1.000 


4.317 


0.915 


0.873 


0.759 


27.20 


19 


5 


0.999 


4.412 


0.970 


0.880 


0.960 


24.52 


21 


4 


1.000 


4.300 


0.966 


0.876 


0.752 


25.75 


29 


4 


1.000 


4.795 


1.160 


0.886 


0.565 


23.4 


29 


5 


0.999 


4.683 


1.132 


0.893 


0.699 


22.5 


29 


6 


1.000 


4.772 


1.164 


0.897 


0.757 


21.9 


29 


7 


0.999 


4.706 


1.160 


0.893 


1.079 


20.9 


39 


4 


0.998 


4.562 


1.006 


0.884 


0.660 


23.6 


39 


5 


1.000 


4.678 


1.026 


0.892 


0.729 


22.3 


39 


6 


1.000 


4.656 


1.020 


0.894 


0.833 


21.8 



TABLE VII. Same as Table but for A 2 = 100/i 2 



/'i 



20/x" 



/'! 



40/i 2 , and \i\ = 60/x 2 



K 


N ± 


(:^(0):> 


( n B,<r) 


(™B,- CT ) 


(n B ) 


(VB,cr) 


{VB,-a) 


(vb) 


11 


4 


0.989 


- 


- 


0.245 


- 


- 


0.1524 


11 


5 


0.998 


- 


- 


0.227 


- 


- 


0.1310 


11 


6 


0.904 


- 


- 


0.206 


- 


- 


0.1207 


13 


4 


1.002 


- 


- 


0.228 


- 


- 


0.1297 


13 


5 


1.000 


- 


- 


0.230 


- 


- 


0.1370 


15 


4 


1.000 


- 


- 


0.205 


- 


- 


0.1149 


15 


5 


1.000 


- 


- 


0.204 


- 


- 


0.1140 


17 


4 


0.998 


- 


- 


0.198 


- 


- 


0.1078 


17 


5 


1.000 


- 


- 


0.212 


- 


- 


0.1203 


19 


4 


1.000 


- 


- 


0.218 


- 


- 


0.1261 


19 


5 


0.999 


- 


- 


0.218 


- 


- 


0.1271 


21 


4 


1.000 


- 


- 


0.215 


- 


- 


0.1245 


29 


4 


1.000 


0.1051 


0.1152 


0.2203 


0.06028 


0.06653 


0.12680 


29 


5 


1.000 


0.0894 


0.1231 


0.2125 


0.04925 


0.07260 


0.12185 


29 


6 


1.000 


0.0798 


0.1317 


0.2115 


0.04334 


0.07794 


0.12128 


29 


7 


1.000 


0.0787 


0.1350 


0.2137 


0.04357 


0.08055 


0.12412 


39 


4 


1.000 


0.1069 


0.1090 


0.2159 


0.06156 


0.06396 


0.12552 


39 


5 


1.000 


0.0907 


0.1232 


0.2139 


0.05109 


0.07272 


0.12381 


39 


6 


1.000 


0.0836 


0.1297 


0.2133 


0.04611 


0.07694 


0.12305 



13 



For higher resolution (K = 21 and 29 for A 2 = 50/i 2 and K = 29 and 39 for A 2 = 100//) 
we plot the structure function f B in Figs. § and ||. Typical contributions to Jb from one- 
boson and two-boson states are shown in Fig. |j. The two-boson contribution to the structure 
function is further analyzed in terms of its dependence on both longitudinal momentum 
fractions in Fig. |^, where 



/b* (2/1,2/2) = J d 2 q ±1 d 2 q ±2 \4> { o 



2,0,0,0) 



(£ 3 )l : 



(3.9) 



is plotted. We also show typical two-body wave functions in Fig. || the agreement with 
the necessary L z — 1 symmetry in the antiparallel helicity case is an important check of J z 
conservation in the calculation. 
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(c) (d) 

FIG. 2. The boson structure function fg at various numerical resolutions for (a) 

(:0 2 (O):) = 0.1, (b) (:</> 2 (0):) = 0.25, (c) (:<p 2 (0):) = 0.5, and (d) (:</> 2 (0):) = 1.0, with M = p, 
A 2 = 50/i 2 , n1 = 10/x 2 , /u| = 20/x 2 , and /x 2 = 30/x 2 . The solid line is from first-order perturbation 
theory. 

As a check on the logarithmic PV coupling constraint in Eq. ( f2.5| ), we compute the bare 
mass shift 5M 2 when M 2 is much less than /z 2 . Figure ^ shows its behavior as a function 
of the longitudinal resolution for various bare couplings. If the logarithmic constraint was 
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FIG. 3. Same as Fig. | but for A 2 = lOOii 2 , n\ = 20ii 2 , t* 2 , = 40tf 2 , and n\ = 60/A 
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FIG. 4. The (a) one-boson and (b) two-boson contributions to the boson structure function 
fs at various numerical resolutions, with M = /i, (:^ 2 (0):) = 1, A 2 = 50/i 2 , /i 2 = 10/i 2 , /i 2 , = 20/i 2 , 
and /i 2 = 30/x 2 . 
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FIG. 5. The two-boson structure function fg a for K = 29 and N± 
(:(f> 2 (0):) = 1, A 2 = 50// 2 , n\ = lO^ 2 , n\ = 20/z 2 , and /x§ = 30^ 2 . 



7, with M 



/i, 





(a) ^ (b) 

FIG. 6. The one-boson amplitude 4> a s ' , with (a) parallel (s = cr) and (b) antiparallel 
(s = — <j) bare helicity, as a function of longitudinal momentum fraction y and one transverse 
momentum component q x in the q y = plane. The parameter values are K = 29, N± = 7, M = //, 



A^ 



50^ 2 , M 2 



10^ 2 , M 2 



20^ 2 , a*§ 



30^ 2 , and (:^ 2 (0):) = 0.25. 
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sufficient for nonperturbative calculation, 5M 2 should go to zero as M 2 approaches zero and 
the resolution approaches the continuum limit. This appears to work well only for weak 
coupling. 
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FIG. 7. The bare mass shift 5M 2 as a function of the longitudinal resolution K for various 



bare couplings g, with N ± = 4, M 2 = O.OOl^ 2 , A 2 = 50/z 2 , fi 2 



W/j 2 , tf. 



20/x 2 , and fij = 30/A 



IV. CONCLUSIONS AND PROSPECTS FOR THE APPLICATION OF 
DLCQ(3+1) TO GAUGE THEORY 



We have used the discretized light-cone quantization method to successfully solve for the 
mass and light-cone wave functions of a dressed fermionic state in a Yukawa theory in 3 + 1 
space-time dimensions. No a priori constraint on the number of bosonic constituents was 
necessary; however, since the fermion constituent was treated as heavy, states containing 
fermion-anti-fermion pairs were truncated. We have found that the eigensolution at strong 
coupling displays features which significantly deviated from first-order perturbation theory. 
Numerical resolution in this domain of strong coupling was not a limiting factor in this 
analysis. 

The regularization procedure which we have chosen, with three PV scalars, functioned 
well. However, better convergence of the DLCQ method at strong coupling could be ob- 
tained by constraining the PV couplings non-perturbatively, rather than using the one-loop 
perturbative constraints in (|2.5|). 

A number of properties of the Yukawa-theory eigensolution could be extracted from 
its light-cone Fock-state wave function. This illustrates the power of the DLCQ method 
in making the Fock-state wave functions of the eigenstates explicitly available. It is also 
possible to use the derived wave functions to compute the Pauli and Dirac form factors of 
the dressed fermion state at general momentum transfers | |29|| . 

There are additional calculations which might be done within the context of the zero 
fermion pair approximation. We could consider two-fermion states and study true bound 
states and scattering solutions. We could also consider dressed spin-3/2 states and the 
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analog of Ntt <-> A transitions. Extension of these methods to pseudo- scalar Yukawa theory 
would make the N7r <-> A connection stronger. 

Although the approach used in this paper to solve (3+l)-dimensional Yukawa theory at 
strong coupling has been successful, future progress for solving other quantum field theo- 
ries will require more efficient analytic methods and numerical algorithms. An alternative 
ultraviolet regularization procedure, using only one PV scalar and one PV fermion |]30|j28[| 
may potentially provide a more efficient approach for solving Yukawa theories. Not only is 
the constraint on the couplings trivial, but the light-cone Hamiltonian is also much simpler. 
The simplifications occur because the instantaneous fermion interactions (the terms in (\Z.'2\) 
of order g 2 ) cancel. Moreover, the DLCQ matrix for the remaining three-point interactions 
is much more sparse, allowing calculations at higher resolutions. Our next step will be to 
test this alternative regularization. 



One can consider two other possibilities for PV regularization |30| of full Yukawa theory. 
One is to use two heavy fermions and one heavy scalar |£|; the other is to use one less 
heavy fermion but make the transverse momentum cutoff part of the regularization rather 
than just a numerical procedure. In each case a 4 term must be added. We plan to explore 
both of these methods. 

Quantum electrodynamics and quantum chromodynamics in physical space-time, in- 
cluding the phenomenon of chiral symmetry breaking remain the central challenge to DLCQ 
methods. One attractive possible approach is to use broken supersymmetry as an effective 
ultraviolet regulator of the light-cone Hamiltonians of gauge theories. 

The PV method also has applicability to the renormalization of non-Abelian Hamiltonian 
gauge theories on the light-cone. Paston, Franke and Prokhvatilov |31].|32j have recently 
extended their analysis to the nonperturbative regulation of light-cone QCD, including the 
regularization of the infrared singularities introduced by using light-cone gauge. They find 
that a combination of light-cone gauge, PV fields, higher derivative regulation, and carefully 
chosen momentum cutoffs can regulate the theory in such a way as to provide agreement with 
Feynman calculations using the Mandelstam-Leibbrandt prescription |33j for the spurious 
singularity in the gauge propagator. The resulting dynamical operator is rather complex 
and several regulating fields are needed, but the regulation procedures appear suitable for 
numerical calculations. We plan to test these methods in Abelian theory, including the 



calculation of positronium bound states 0,0 and the non-perturbative calculation of the 
anomalous magnetic moments of leptons at strong coupling |36] . 
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APPENDIX: LANCZOS ALGORITHM FOR INDEFINITE METRIC 



The ordinary Lanczos algorithm 27] was designed for diagonalization of real symmetric 



or Hermitian matrices. A more general form, the biorthogonal Lanczos algorithm |3"7| , can 



be applied to non-symmetric cases but is quite cumbersome. In the case of a complex 
symmetric matrix, the biorthogonal algorithm can be reduced to a form nearly as simple as 
the real symmetric case [BS| ; this approach was used in previous work fl8| , [5If where imaginary 
couplings made negative norms unnecessary. The complex-symmetric approach is not easy to 
implement for Yukawa theory because the Hamiltonian is fully Hermitian. Instead, negative 
norms are assigned, and the eigenvalue problem becomes one with indefinite metric. 

For this case the biorthogonal algorithm can still be reduced to a simpler form. Let rj 
represent the metric signature, so that numerical dot products are written {<p'\<p) = </>'* ■fjtp- 
The Hamiltonian matrix H is constructed to be self-adjoint with respect to this metric, 



which means that [^] H 
H then takes the form 



H. The Lanczos algorithm for the diagonalization of 



Q i 



lj+l 



v 3 q] ■ riEqj , Tj = Hq j - Ti-i^-i 



Va 



.1 = signfr • ■ T]rj 



3) ' 



v\ 



- ^jQj , Pj = 
sign(qr* ■ rjq^ 



lj 



'iFi\ 



■j+lVjPj , 



(Al) 



v 



where q l is taken as a normalized initial guess and 70 = 0. This initial guess is generated 
with use of high-order Brillouin-Wigner perturbation theory. To determine when to stop 
the Lanczos iterations, the convergence of the eigenvalue and parts of the boson-fermion 
wave function are monitored. 

Just as for the ordinary Lanczos algorithm, the original matrix H acquires the following 
tridiagonal matrix representation T with respect to the basis formed by the vectors q ,: 



H -^T 



/«i /3i 

7l «2 /?2 
72 "3 P3 

73 • ■ 
0.. 

V 



(A2) 



By construction, the elements of T are real. The new matrix is also self-adjoint, but with 
respect to an induced metric v = {vi,V2, ■ ■ ■}■ The eigenvalues of T approximate some of 
the eigenvalues of H, even after only a few iterations. Approximate eigenvectors of H are 
constructed from the right eigenvectors q of T as 4> i = J2k( c i)kQk- 
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